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We present an update of the measurement of Qs(m^) from ALEPH r hadronic spectral func- 
tions. We report a study of the perturbative prediction(s) showing that the fixed-order per- 
turbation theory manifests convergence problems not presented in the contour-improved cal- 
culation. Potential systematic effects from quark-hadron duality violations are estimated to 
be within the quoted systematic errors. The fit result is as{m%) = 0.344 ± 0.005 ± 0.007, 
where the first error is experimental and the second theoretical. After evolution, the as{m%) 
determined from r data is the most precise one to date, in agreement with the corresponding 
N^LO value derived from Z decays. 



1 Introduction 

The T lepton, through its hadronic decays, provides a clean laboratory to perform precise studies 
of QCD. Invariant mass distributions obtained from long distance hadron data allow one to 
compute the spectral functions, which permit the study of short distance quark interactions. In 
particular, these spectral functions can be exploited to precisely determine the strong coupling 
constant at the r-mass scale, as(m^). The present analysis is described in detail in re 

2 Tau Hadronic Data and Spectral Functions 

The nonstrange vector (axial- vector) spectral functions vi{ai), for a spin 1 hadronic system, are 

obtained from the squared hadronic mass distribution, normalised to the hadronic branching 
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fraction (with Rry/A = ^^''"^ )> ^^'^ divided by a factor exhibiting kinematics and spin 

characteristics 
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The basis for comparing a theoretical description of strong interaction with hadronic data is 
provided by the optical theorem, which relates the imaginary part of the polarisation functions 
on the branch cut along the real axis, to the spectral functions: 27r • Imlly^^{s) = vi/ai{s). 

The total hadronic observable Rr is obtained from measured leptonic branching ratios, or 
only from the electronic one assuming universality. The two determinations are in very good 
agreement, yielding Rr = {I - Be - B^)/Be = l/Bf' - 1.9726 = 3.640 ±0.010. One can identify 
in Rr a component with net strangeness and two nonstrange vector(V) and axial-vector(A) com- 
ponents. Including the latest results from BABAR and Belle the value of the strange component 
is Rr,s = 0.1615 lb 0.0040 . The separation of the V and A components is straightforward for 
final states with only pions using G-parity. However, KK modes are generally not eigenstates 
of G-parity. The decay to K~K^ is pure vector. The vector component of the KKn mode 
is determined assuming CVC and using new measurements from the BABAR Collaboration^, 
for the e~^e~ annihilation to K^K~tt^ and to K^K^tt^. After integration one gets a clear 
dominance of axial- vector component, fA,CYc{KKT:) = 0.833 it 0.024. For the KKtttt rarer 
modes a conservative value fAiKK-Kir) = 0.5 ± 0.5 is used. Finally, we get the components: 
Rry = 1-783 ± 0.011 ± 0.002 and Rr,A = 1-695 ± 0.011 ± 0.002 , where the first errors are 
experimental and the second due to the V/A separation. 

3 Theoretical Prediction of Rr 

The nonstrange ratio Rry /a can be written as an integral of the spectral functions over the 
invariant mass-squared s of the final state hadrons 
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The two point correlator can not be predicted by QCD in this region of the real axis. However, 
using Cauchy's theorem, one can relate this expression to an integral on a circle in the complex 
plane. Then, the OPE yields 



Rry/A cx 1 + + 5W + ^Sy/l + T.€y/A ' (3) 

D=4,6,, 



with a massless perturbative contribution, a non-logarithmic electroweak correction, the dimen- 
sion two perturbative quark-mass contribution and higher dimension nonperturbative conden- 
sates contributions respectively. The perturbative part reads 6^'^^ = J2'?^=i KniO^^'^H'^s) , with 
the functions 

^<..K, = ^ ft(r-2±,2(±]'-(±]'](^y, (4) 
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where ^ is a scale factor. A breakthrough was made recently^, so that the pertubative coefficients 
are now known up to K4 (see^ for the numerical values of the Kn {0 coefficients) . 



3.1 Perturbative Methods 

The perturbative contribution to Rr provides the main source of sensitivity to as(so). The 
value of the strong coupling in the complex plane can be computed assuming the validity of 
the renormalisation group equation (RGE) outside the real axis, and using a Taylor series of 
r] = ln{s/so). In the fixed order perturbation theory (FOPT), at each integration step, the Taylor 
expansion is made around the physical value as(so)- This may cause important problems as the 



absolute value of rj gets large and the convergence speed of the series is reduced^ . In addition, 
a cut at a fixed order in as(so) is applied on the Taylor series and on the integration result 
in FOPT. Therefore, important known higher order terms are neglected, yielding additional 
systematic uncertainties. A better suited method is CIPT which, at each integration step, 
computes as{s) using the value found at the previous step. In this approach the Taylor expansion 
is always used for small absolute values of its parameter, hence excellent convergence properties. 

In practice we have also used geometric growth estimations for the first unknown coefficients 
Pi, and K^. We have tested that CIPT is less sensitive to changes of these coefficients, 
and it also exhibits a smaller scale dependence than FOPT. Numerically, the difference of the 
perturbative contributions computed with the two methods are about 15%. In fact this difference 
could have been much larger if not for the properties of the kernel in the integral (jH) which has 
small absolute values in the region where the as{s) predictions of the two methods are rather 
differentia. 

The CIPT method behaves better than FOPT and is to be preferred. The difference between 
the results obtained with the two approaches is not to be interpreted as a systematic theoretical 
error, but rather like a problem of foptE. 

3.2 Quark-Hadron Duality Violation 

It is known that OPE describes only part of the nonperturbative effects. In order to estimate the 
impact of the missing contributions, we test two models based on resonances and on instantons. 
We add their contributions to the theoretical prediction, choosing parameters that provide a 
good matching to the V+A spectral function near the r mass. For these models, we find 
corrections situated within our systematic uncertainties^. 

4 Combined Fit 

In order to obtain additional experimental information, we use spectral moments defined as 



They allow one to better exploit the shape of the spectral functions and they suppress the 
region where OPE fails. The corresponding theoretical prediction is very similar to ([3]), with 
consequent perturbative and nonperturbative contributions. Due to strong correlations, we 
use only Rr {k = and I = 0) and four additional moments (A: = 1 and / = 0,1,2,3) to 
simultaneously fit as(m^) and the leading D = 4, 6, 8 nonperturbative contributions. 

In spite of the fact that the nonperturbative contributions fitted for the V and A spectral 
functions have opposite signs and they are one order of magnitude larger than those from V+A, 
we find an excellent agreement between the values found for ag(m^) from the three fits. The 
result of the fit to the V+A spectral moments reads 



where the first error is experimental and the second is theoretical. When evolving this value to 
the Z scale^I^see Fig. [TJ one gets 





as{ml) = 0.344 ± 0.005 ± 0.007 



(6) 



a^J'>{ml) = 0.1212 ± 0.0005 ± 0.0008 ± 0.0005 



(7) 



where the first two errors are propagated from ([6]), and the last one summarises uncertainties 
in the evolution. The consistency between this result and the value found by a global fit to 
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Figure 1: Top: The evolution of asim^) to higher scales fi using the four-loop RGE and the three-loop matching 
conditions applied at the heavy quark-pair thresholds (hence the discontinuities at 2mc and 2m(,). The evolution 
is compared with independent measurements covering /x scales that vary over more than two orders magnitude. 
Bottom: The corresponding as values evolved to mz- The shaded band displays the r decay result within errors. 

electroweak data at the Z-mass scale^, a^g\m^) — a^^^m^) = 0.0021 it 0.0029, provides the 
most powerful present test of the evolution of the strong interaction coupling over a range of s 
spanning more than three orders of magnitude. 

5 Conclusions 

Motivated by some new results both on theoretical and experimental grounds, we have revisited 
the determination of as(m^) from the ALEPH r spectral functions. We have reexamined two 
common numerical methods: we have identified specific consistency problems of FOPT, which 
do not exist in CIPT. The r measurement of evolved to the Z scale is found to be in excellent 
agreement with the direct determination from Z decays. Both results are the only ones at N^LO 
order so far, confirming the running of as between 1.8 and 91 GeV, as predicted by QCD, with 
an unprecedented precision of 2.4%. 
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